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ABSTRACT 

Weakly interacting massive particles (WIMPs) may be indirectly detected by observation of upward muons 
induced by energetic neutrinos from annihilation of WIMPs that have accumulated in the Sun and/or Earth. 
Energetic muon neutrinos come from the decays of r leptons, c, 6, and t quarks, gauge bosons, and Higgs bosons 
produced by WIMP annihilation. We provide analytic expressions, suitable for computing the flux of upward 
muons, for the neutrino energy spectra from decays of all these particles in the center of the Sun and Earth. 
These analytic expressions should obviate the need for Monte Carlo calculations of the upward-muon flux. We 
investigate the effects of polarization of the gauge bosons on the neutrino spectra, and flnd that they are small. 
We also present simple expressions for the second moments of the neutrino distributions which can be used to 
estimate the rates for observation of neutrino-induced muons from WIMP annihilation. 
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1. Introduction 

Weakly-interacting massive particles (WIMPs) with masses in the GeV- 
TeV range are among the leading candidates for the dark matter in our Galactic 
halo [1]. One of the most promising methods for discovery of WIMPs in the 
halo is via observation of energetic neutrinos from annihilation of WIMPs in 
the Sun and/or Earth [2] [3] [4] [5]. WIMPs in the halo can be captured in the 
Sun and the Earth by elastic scattering from nuclei therein. WIMPs that have 
accumulated in the Sun or Earth will annihilate into ordinary particles, such 
as quarks, leptons, and if heavy enough, gauge bosons, Higgs bosons, and top 
quarks. Many of the leading WIMP candidates — such as the neutralino in 
supersymmetric theories — are Majorana particles, so they do not annihilate 
directly into neutrinos [6]. Neutrinos with energies in the multi-GeV range are 
then produced by decay of the annihilation products, and such neutrinos are 
potentially observable in many astrophysical neutrino detectors. 

The differential energy flux of high-energy neutrinos of type i (e.g., i = 
VfijVfi, etc.) from WIMP annihilation in the Sun (Earth) is 
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where F^ is the rate of WIMP annihilation in the Sun (Earth), and R is the 
distance of the Earth from the Sun (or the radius of the Earth). The sum is 
over all annihilation channels F (e.g., pairs of gauge or Higgs bosons or fermion- 
antifermion pairs), Bp is the annihilation branch for channel F, and (dN/dE)Fi 
is the differential energy spectrum of neutrino type i at the surface of the Sun 
expected from injection of the particles in channel F in the core of the Sun. The 
total number of neutrinos of type i from final state F is Np^i = f(dN/dE)pidE. 
The spectrum (dN/dE)Fi is a function of the energy of the neutrino and of the 
energy of the injected particles. As we point out in this paper, the (dN/dE)Fi 
may also depend on the polarization of the annihilation products. 

This paper continues the work of Gaisser, Steigman, and Tilav [4] and 
Ritz and Seckel [3]. The dN/dE from the r-lepton and h- and c-quark final 



states were first computed analytically using the spectator-quark approxima- 
tion in Ref. [4], but stopping of heavy hadrons and stopping and absorption 
of neutrinos in the Sun were neglected. A Monte Carlo calculation which in- 
cluded interactions with the solar medium was performed by Ritz and Seckel 
[3] (hereafter RS). They also considered decays of the top quark, but those 
results, obtained assuming a top-quark mass of 60 GeV, are no longer valid. 
Some of the recent experimental collaborations have followed the example of 
RS and used Monte Carlo simulations to obtain the neutrino energy spectra 
[7] [8]. However, Monte Carlo calculations can be cumbersome, especially when 
numerous WIMP candidates are to be explored; analytic expressions are much 
more convenient. It is also easier to isolate the various interaction effects and 
theoretical assumptions that enter the calculation with analytic expressions. In 
this paper, we derive analytic expressions for the neutrino energy distributions 
from r-lepton, c-, &-, and t-quark, and gauge- and Higgs-boson decay in the Sun 
and Earth that can be used in place of Monte Carlo results. 

Our results for h and c quarks and r leptons are derived from the results 
of [3] and [4]. Our results for dN/dE from gauge-boson and Higgs-boson final 
states follow the work of Refs. [9] and [5], but our discussion of the top-quark 
final state is entirely new. If the top quark is heavier than the W boson, 
then it decays almost exclusively into W bosons and h quarks; therefore, the 
calculation of the neutrino spectrum must be redone. Also note that in many 
cases where the neutralino is heavier than the top quark, neutralinos annihilate 
predominantly to the tt final state [10] [11] [12]. Therefore, evaluation of the 
neutrino spectrum for the ti final state is crucial to calculation of energetic- 
neutrino event rates. 

Of all our results, the neutrino spectra from injected h and c quarks are 
subject to the largest theoretical uncertainty. In particular, heavy hadrons 
will be slowed before they decay as they pass through the Sun [3], and little 
is known about interactions of heavy hadrons with a surrounding medium. 
However, this uncertainty enters the Monte Carlo calculations as well, so little 
accuracy is lost by using analytic results instead. In addition, the effect of 
heavy-quark stopping becomes small at low injection energies, so the resulting 
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uncertainty is small at low energies. Once the neutralino mass has become large 
enough that the top-quark annihilation channel opens up, the branching ratio 
for neutralino annihilation into h quarks becomes negligible. So although the 
&-quark neutrino spectrum may be quite uncertain at high injection energies, 
the resulting uncertainty in the total neutrino event-rate remains small. Note 
that in all cases, the errors in the upward-muon flux from the neutrino spectra 
are dwarfed by underlying uncertainty in the local halo density. 

When the WIMP is heavy enough that it annihilates primarily to vector 
bosons and/or top quarks, then the theoretical calculation simplifles greatly. 
These particles (as well as r leptons) do not hadronize and are not slowed 
before they decay, so relatively simple analytic expressions are obtained for the 
neutrino spectra from the Earth. The only point at which the analysis becomes 
somewhat complicated is when considering the stopping of the h quark from 
top quark decay in the Sun. However, this &-quark contribution provides a 
small fraction of the total neutrino signal from top decay, so the uncertainty 
introduced from modeling of the hadronic stopping effect is negligible. In the 
end, we arrive at expressions for dN/dE from the Earth which give results for 
the upward-muon flux estimated to be accurate to a few percent and results for 
dN/dE from the Sun which give results accurate to within 5%. 

We also calculate the neutrino spectra from polarized gauge bosons. In Ref. 
[5] it was assumed that the gauge bosons produced by neutralino annihilation 
were unpolarized. In fact, they are indeed polarized [13]. In this work, we 
quantify the effect of polarization of gauge bosons from neutralino annihilation 
on the upward-muon flux from decays of these particles and show that it leads 
to no more than a 10% correction to the results of Ref. [5]. 

In most neutrino detectors, the existence of the neutrino is inferred by ob- 
servation of an upward muon produced by a charged-current interaction in the 
rock below the detector. The cross section for production of a muon is pro- 
portional to the neutrino energy, and the range of the muon in rock is roughly 
proportional to the muon energy. Therefore, the rate for observation of ener- 
getic neutrinos is essentially proportional to the second moment of the neutrino 
energy spectrum. Inserting the numerical values for the charged-current cross 



section and the effective range of the muon, and ignoring detector thresholds, 
the rate per unit detector area for neutrino-induced throughgoing-muon events 
may be written [3] [5] 



rdetect = (1.27 X 10-29 m-2yr-l) 
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(1.2) 
for neutrinos from the Sun; the same expression multiplied by 5.6 x 10® (the 
square of the ratio of the Earth-Sun distance to the Earth's radius) gives the 
rate for neutrino events from the Earth. Here, C is the WIMP capture rate in 
the Sun or Earth, Bp is the branching fraction for WIMP annihilation to the 
flnal state F, the sum on i is over muon neutrinos and muon anti-neutrinos, 
the flj- are neutrino scattering coefficients, a^ = 6.8 and ap = 3.1, and the hi are 
muon range coefficients, hi, = 0.51 and hp = 0.67. The quantity 
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is the second moment of the spectrum of neutrino type i from final state F 
scaled by the square of the injection energy Ei of the annihilation products. 
The quantity z = E,^/Ei„ is the neutrino energy scaled by the injection energy. 

Strictly speaking, neutrino telescopes observe neutrinos only with ener- 
gies above a given threshold. Therefore, the event rate is proportional to 
the contribution to the second moment from neutrinos with energies above 
threshold — that is, there is a lower bound to the integral in Eq. (1.3). To 
obtain the most accurate experimental information, a detailed calculation of 
the neutrino spectra should be folded in with the detector response. For ex- 
ample, an energetic-neutrino signal from the Sun must be distinguishable from 
the atmospheric-neutrino background. The direction of the neutrino-induced 
muon is correlated with the parent-neutrino direction only within an angular 
window of roughly ^jGeW / Ei, radians. Therefore, a proper determination of 
the signal-to-noise ratio requires knowledge of the neutrino energy spectrum. 

On the other hand, in cases where the WIMP is quite massive, most of 
the neutrinos have energies large enough to produce a muon with an energy 



much higher than the detector threshold. For example, if a 100-GeV WIMP 
annihilates to gauge bosons, the typical neutrino energy is roughly half that, and 
the typical muon energy typically half that, about 25 GeV. Typical thresholds 
for current and next-generation detectors are no more than 10 GeV, so the vast 
majority of the neutrinos in this example are above threshold. In addition, even 
if a non-negligible fraction of the neutrinos have energies below the detector 
threshold, the total second moment of the energy distribution is still primarily 
determined by the higher-energy neutrinos. The contribution of neutrinos with 
energies below threshold to the upward- muon flux will be C'(-E'fhresh/"^x^) where 
-E'thresh IS the threshold energy. Therefore, in cases where the WIMP mass 
is significantly larger than the detector threshold, the expression for the rate 
for neutrino-induced upward muons, Eq. (1.2), together with the results for 
i^N z'^j presented here will provide a good theoretical estimate to compare with 
experimental determinations of the fiux of upward muons from the Sun and/or 
Earth. 

The most promising method of detection of the energetic neutrinos, es- 
pecially for higher-mass WIMPs, is the upward-muon signal. Therefore, we 
calculate the energy spectrum of muon neutrinos and antineutrinos only. 

In the next Section, we discuss some general results on neutrino spectra. 
We relate the neutrino spectrum from decay of a moving particle to the neutrino 
spectrum from decay of the particle at rest, and we show how the neutrino 
spectra from particle decays in the Sun are obtained from those from decay 
in the Earth. We discuss how hadronization and stopping of h and c quarks 
affect the neutrino spectra. We also discuss how neutrino spectra from Higgs 
bosons are obtained from the neutrino spectra from the Higgs-boson decay 
products. Much of this is review, but useful for the discussion in the Sections 
that follow. In Section 3, we calculate the neutrino spectra from decays of h 
and c quarks and r leptons. In Section 4, we discuss the neutrino spectra from 
decays of gauge bosons and quantify the effect of gauge-boson polarization on 
the neutrino spectra. In Section 5, we calculate the neutrino spectra from 
decays of top quarks, and we summarize our results in Section 6. 



2. Preliminaries 

For all of the particles considered, calculation of the neutrino distribution 
in the rest frame of the decaying particle, {dN / dEf^^^ {E u) , as a function of the 
neutrino energy E^, is a standard exercise. Once the rest-frame distribution is 
known, the energy distribution of a particle moving with an energy Ed, velocity 
/3, and 7 = (1 — /3^)~"'"" = Ed/rrid (where rrid is the decaying-particle mass) is 
related to the rest-frame distribution by 
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This assumes that there is no prefered direction in the rest frame of the parent 
particle. In particular, it applies to decay of scalar or unpolarized particles. If 
the decaying particle is polarized and the neutrino direction is correlated with 
the spin of the decaying particle, Eq. (2.1) does not apply. In almost all cases 
we consider, the decaying particles are unpolarized, and we show that in the 
few cases where they are polarized, the effect of polarization is small or can be 
accounted for relatively easily. 

The Earth is thin enough that stopping of heavy hadrons and stopping 
and absorption of neutrinos as they pass through the Earth can be neglected. 
Therefore, if the rest-frame distribution from a particular decaying particle is 
known, Eq. (2.1) gives the correct dN/dE for particle decay in the Earth; 
hence the superscript © in Eq. (2.1). This also implies that the neutrino and 
antineutrino spectra from injection of a given particle-antiparticle pair in the 
center of the Sun are the same. Also, if the scaled second moment in the rest 
frame of the decaying particle, i^N z'^j , is known, then the scaled second 
moment for a particle that decays with a velocity (3 in the Earth is simply 



(#z2)(/3) = (#z2p(l + ^). 



(2.3) 
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On the other hand, energetic neutrinos will lose energy via neutral-current 
interactions with the solar medium and become absorbed via charged-current 
interactions as they pass through the Sun. To account for these effects, we use 
the calculation of RS that a neutrino injected with an energy E leaves the Sun 
with energy 

(2.4) 
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where r^ = 1.01 x IQ-^ GeV"! and rp = 3.8 x lO"-* GeV-\ and with probability 

1 



P 



I 



1 + En 



(2.5) 



where a^ = 5.1 and ap = 9.0 As a result, the neutrino spectrum for a particle 
decaying with energy Ed in the Sun, (dN/dE)'^ , is related to the neutrino 
spectrum for a particle decaying with energy Ed in the Earth, (dN/dE)® , by 
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where E^. = E^/(l — E^n) is the energy a neutrino had at the core of the Sun 
if it exits with energy Ei, . Note that stopping and absorption are different for 
neutrinos and antineutrinos, so the spectrum of neutrinos from WIMP annihi- 
lation in the Sun is different than the spectrum of antineutrinos. The spectrum 
of neutrinos is the same as that for antineutrinos for decays in the Earth, so we 
neglect the subscript i on neutrino spectra from the Earth. 

If a & or c quark is injected into the center of the Earth, it will lose energy 
during hadronization. As a result, the energy at which the hadron decays in the 
Earth, Ed, is related to the energy Ei at which it is injected by Ed = zjEi. We 
are interested primarily in the second moment of the neutrino distributions, so 
we use the square root of the second moments of the fragmentation functions 
computed by RS — that is, we take zj = 0.58 for c quarks and zj = 0.73 for h 
quarks. Actually, the distribution oizj is described by a fragmentation function 
[14], but this fragmentation function is highly localized around zj . Comparison 
of higher order moments of zj suggests that no more than 0(5%) accuracy is 
lost by using the central value of the fragmentation function. 



In addition, the core of the Sun (but not the Earth) is dense enough that 
h- and c-quark hadrons will interact with the solar medium and be slowed 
appreciably before they decay [3]. If a hadron initially has an energy Eq = zjEi 
(after hadronization of a quark injected with energy Ei), then it will decay with 
an energy Ed picked from a decay distribution [3], 
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where Eg = 250 GeV for c quarks. Eg = 470 GeV for h quarks. 

The nth moment of this distribution, for a given initial energy Eg, is 
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so in particular, the average energy at which a hadron decays is 



{Ed)iEo) = Er,exp[^]EJ^^ 



(2.8) 



(2.9) 



where Ei(x) = J (e ^ ly)dy [15], and the rms value of the decay energy is 

ET%Eo) = ^) = VEciEo-{Ed)). (2.10) 

These quantities will come in handy in the following Section. 

Of all the factors that enter into the calculation of the neutrino spec- 
trum, stopping of heavy hadrons is the most uncertain theoretically; very lit- 
tle is known about the interaction of heavy-quark hadrons with dense matter. 
The functional form of the decay distribution, Eq. (2.7), is physically well- 
motivated, but the values of Eg listed here are subject to some (perhaps size- 
able) uncertainty. However, this theoretical uncertainty also enters the Monte 
Carlo calculation; in this regard, Monte Carlo simulations offer no improvement 
over the analytic result. Also, the effect of hadron stopping is small at lower 
energies, so the uncertainty introduced into the neutrino spectra is relatively 
small. Stopping becomes much stronger — and the subsequent theoretical un- 
certainty much larger — at higher injection energies, but in most cases where 
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the WIMP is massive enough that it can annihilate into top quarks, it anni- 
hilates almost exclusively to top quarks, gauge bosons, and/or Higgs bosons 
[10][11][12]. Therefore, the uncertainty in the total neutrino spectrum due to 
poor understanding of hadron stopping is never very large. 

Before continuing, we remind the reader that light (u, d, and s) quarks will 
produce long-lived hadrons that come to rest in the core of the Sun or Earth 
before they decay; therefore, these final states produce no energetic neutrinos. 
The same is true of neutrinos from muon decay. On the other hand, r leptons, 
top quarks, and gauge and Higgs bosons will decay essentially immediately, so 
there is no energy loss due to slowing or hadronization before they decay. 

Finally, consider a chain decay process, H ^ I1I2, where a heavy particle, 
H, decays to two lighter particles, /i and I2, which then decay to energetic 
neutrinos. For example, the heavy particle could be one of the Higgs bosons 
which then decays to two light fermions. In the rest frame of the decaying H , 
the energies of the two light particles are Ei = {mjj -\- m\ — m^) / {2mH) and 
E2 = {rn^H + m\ — m\)/(2mH), and their velocities are /3i = (1 — m'f/E'fy'''^ 
and /32 = (1 — m'^/E'^YI'^, where niH , rrii, and 1112 are the masses oi H , li, and 
I2, respectively. Suppose the neutrino spectrum from decay of /i as a function 
of neutrino energy Ei, for a given injection energy Ei is (dN/dE)i(Ei, E^), 
and similarly for neutrinos from decays of /2- Then the spectrum from the 
heavy particle H, which decays while moving with a velocity i3h and energy 
Eh = InmH is 
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where the sum is over the two decay particles, and Ej and (3j are the decay- 
particle energies and velocities in the rest frame of the decaying H . Also, 
if the scaled second moments of the neutrino distribution, i^N z'^j {Ei) and 
(Nz'^\{Ei), from decays of /i and I2, respectively, are known for injection 
energies Ei, then the i^N z'^j from decay oi H ^ I1I2 is [16] 
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These equations can be used, for example, to determine the neutrino spectra 
from Higgs-boson decay. Higgs bosons do not decay directly to neutrinos, so 
all the energetic-neutrino signal comes from the decay of the Higgs decay prod- 
ucts. Given the neutrino spectra from the Higgs decay channels (e.g., quarks, 
leptons, and if heavy enough, gauge bosons and/or top quarks) presented be- 
low, the neutrino spectrum and/or i^N z"^) from Higgs decay are obtained by 
summing Eqs. (2.11) and/or (2.12) over all decay channels of the Higgs with 
the appropriate branching ratios. We will use these results below to obtain the 
i^N z"^) from top-quark decay. Also, note that these relations, Eqs. (2.11) and 
(2.12), are valid for neutrino spectra from particle decay in both the Sun and 
the Earth. 

Given the results in this Section, the neutrino spectra from each final state 
follow from the kinematics of two- and three-body decay, and we will discuss 
each final state below. 



3. Neutrinos From r Leptons and h and c Quarks 

3.1. Neutrinos From r Leptons 

We begin with the neutrino energy spectra from r leptons. Calculation 
of these spectra is easier than that for h and c quarks, since r leptons do not 
hadronize, and they are not slowed before stopping. In addition, we consider 
decay of unpolarized r leptons only. For neutralinos, this is valid because 
the r leptons produced by neutralino annihilation are unpolarized [13]. If for 
some unforeseen reason, the injected r leptons are polarized, the results for the 
neutrino spectra will be slightly different. 

In the rest frame of the r, the energy distribution of the muon neutrino 
from the decay r -^ jJ-Vf^VT is [17] 
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where Tr^fi^v — 0.18 is the branching ratio for r decay to muons. The spectrum 
of neutrinos from decays of r leptons moving with an energy Ei (the injection 
energy), velocity /3, and 7 = (1 — /3^')~^'^ is 
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where E± = £"^7 "'^(1 =F /?) ^ ■ In the relativistic limit (/3 -^ 1), this becomes 
[dN/dE) = (2r/£'i)(l - 3a;2 + 2a;3)^ for < E"^ < £",•, where x = EjEi. Given 
the neutrino distribution from r decay in the Earth, the spectrum from r decay 
in the Sun is obtained by using Eq. (2.6) to account for stopping and absorption 
of neutrinos. 

Applying the results of Section 2, the scaled second moment from injection 
of TT pairs with velocity (3 in the core of the Earth is 
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The integral over the distributions from the Sun can also be performed, and 
the result is 
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where the r,- are the neutrino stopping coefficients given in Section 2. For muon 
neutrinos, the function hji{y) is given by 
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and for antineutrinos, the appropriate function is 
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In deriving Eq. (3.5), we took a^ = 5.0, which is an excellent approximation to 
the RS value of a^ = 5.1. In addition, the h functions were obtained using the 
j3 ^t \ Winiioi {dN / dE)f^. If /3 is not near unity, then Ei must be close to m^, in 
which case neutrino absorption and stopping will be negligible and the dN/dE 
from the Earth can be used. Also, in the low-energy limit {Ei ^ r~ ~ TeV), 



simple but accurate expressions can be obtained from Taylor expansions of the 
h functions. 

In addition to the energetic muons produced by direct decays of the r, there 
will be additional indirect energetic muon neutrinos from the r neutrinos. The 
r neutrino from the r decay may produce an additional r lepton via a charged- 
current interaction with the solar medium as it passes through the Sun, and 
this r lepton will then decay and produce an additional muon neutrino [3]. The 
energy of this muon neutrino is much smaller than that from direct decay of the 
original r lepton. Although these neutrinos may contribute non-negligibly to 
dN/dE at low neutrino energies, this contribution to the upward-muon signal 
is tiny, so it is justifiably ignored when considering upward muons. Therefore, 
Eqs. (3.3), (3.4), (3.5), and (3.6) account for the overwhelming majority of the 
energetic neutrinos that give rise to upward muons. In addition, they correctly 
take into account all the important physical processes: the r-decay kinematics, 
and stopping and absorption of neutrinos as they pass through the Sun. There- 
fore, these expressions are subject to little theoretical uncertainty and can be 
used in place of the results of Monte Carlo calculations. 

3.2. Neutrinos from h and c Quarks 

The three-body decay of h- and c-quark hadrons is similar to that for r 
decay, and the treatment of neutrino stopping and absorption in the Sun is sim- 
ilar; however, calculation of the neutrino spectra is slightly more complicated 
due to hadronization of the quarks and stopping of the heavy hadrons in the 
Sun. Still, we can treat these effects analytically without resorting to Monte 
Carlos. Hadronization is simple: the energy of a quark is degraded from its in- 
jection energy Ei to the hadron energy Eg = zjEi. Stopping is straightforward 
as well: we must integrate the neutrino energy spectrum from a hadron which 
decays with energy Ed over the distribution of decay energies, Eq. (2.7). 

In the rest frame of the &-quark hadron, the energy distribution of the 
neutrino from the decay h -^ cjiv is (neglecting the small correction due to 
finite c-quark mass) similar to that from r-lepton decay [17]: 
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where Fj^^^x — 0.103 is the branching ratio for inclusive semileptonic decay 
of the h quark into muons [18]. A h quark injected with energy Ei decays in the 
Earth with an energy Ed = ZfEi, velocity /3, and 7 = (1 — /3^)~"'"" = Ed/rrn, 
and the resulting neutrino spectrum is 
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where E± = £"^7 "'^(1 =F /?) ^ ■ In the relativistic limit (/3 -^ 1), this becomes 
[dN/dE) = (2r/£'d)(l - 3a;2 + 2x^), for < E"^ < -Ed, where x = E^/Ed- 

The neutrino spectrum from h decay in the Sun is obtained by applying 
Eq. (2.6) to Eq. (3.8), and integrating over the distribution of hadron decay 
energies, Eq. (2.7): 
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+ A8y^x + 2y^x'^ - l^^^y'^x - iy'^x^), 

(3.9) 

where E^ = E^/(l — E^Ti), x = inax.(Ec/ Eg, Eg/ E^), Eg = zjEi, and here, 
y = E^/ Eg. The value of Eg to be used here is 470 GeV. Eq. (3.9) was obtained 
using the /3 ^ 1 limit of {dN / dE)®^^, so it is not strictly valid for Ei near toj. 
On the other hand, if the injection energy is so low that (3 is not near unity, 
interactions with the solar medium are small and the spectrum from decay in 
the Earth can be used. 

In addition to the prompt neutrinos from &-quark decay, there will be 
additional energetic neutrinos from subsequent decays of the c quark. How- 
ever, the typical energy of these neutrinos is roughly 1/3 the energy of the 
neutrinos from prompt decay, so they contribute only about 10% to the total 



upward- muon signal. Since neutralino annihilation into h quarks is always ac- 
companied by annihilation into r leptons, c quarks, and possibly top quarks 
and gauge/Higgs bosons, the error introduced into the upward-muon rate by 
including only prompt neutrinos from h decay will never be greater than 10%. If 
greater accuracy is desired, the contribution from c quarks can be included using 
the c-quark neutrino spectra below and the chain-decay equations of Section 2. 
The neutrino spectra from charmed hadrons are similar, but the three-body 
decay kinematics is slightly different. In the rest frame of the charmed hadron, 
the energy distribution of the neutrino from the decay c -^ sjiv is (neglecting 
the small correction due to finite strange-quark mass) [17] 

rest 
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Q<E,<-m,, (3.10) 



where Fc^^^x — 0.13 is the branching ratio for inclusive semileptonic decay of 
the c quark into muons. (We obtain this value for F^^^^x by assuming charged 
and neutral charmed mesons are produced with equal probability and averaging 
the branching ratios for inclusive semileptonic decays to muons [18]. Our value 
for this branching ratio appears to be almost twice as large as that used by 
RS, so our neutrino spectra from c-quark decay will be significantly larger than 
those obtained by RS.) A c quark injected with energy Ei decays in the Earth 
with an energy Ed = ZfEi, velocity /3, and 7 = (1 — /3^)~"'"" = Ed/rric and the 
resulting neutrino spectrum is 
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where E± = E,jj~^{l =p l3)~^. In the relativistic limit {j3 -^ 1), this becomes 
{dN/dE) = {T/Ed)[{5/3)-3x^ + {i/3)x^], for < E^, < Ed, where x = EjEd- 
The derivation of the neutrino spectrum from c decay in the Sun is similar 
to that for h quarks above. The result is 
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(3.12) 
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where as before, Em = E^/(l — E^Ti), x = inax.(Ec/ Eg, Eg/ E^), Eg = zjEi, 
and here, y = E^/ Eg. The value of Eg to be used here is 250 GeV. 

Now we list expressions for the i^N z'^j . The expression for decay of a & 
quark moving with velocity (3 in the Earth is 
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and that from c-quark decay is 
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(1 + y). 



(3.13) 



(3.14) 



The second moments from h- and c-quark decay in the Sun can be obtained 
by integrating the correct expressions for dN/dE. However, it turns out to 
be easier to switch the order of integration — that is, first compute (^Nz^j for 
a given decay energy, Ed, and then integrate over the hadronic decay-energy 
distribution. The result is 
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where the subscript / denotes & or c quarks, and the subscript i refers to neutrino 
type {Vfi or i>^), and the moments {E^) are those that were given in Section 2. 
The integral in the first line cannot be performed analytically, but we have 
found that the approximation in the second line is accurate to a few percent for 
injection energies less than a few TeV. For h quarks, the functions h},i are the 
same as those for r leptons given in the previous subsection: hi,i(y) = hji{y) 
[c.f., Eqs. (3.5) and (3.6)]. For neutrinos, the c-quark functions are, 
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Eqs. (3.13), (3.14), and (3.15), with the expressions for the appropriate h(y) 
functions, correctly describe the second moments of the neutrino distributions 
from h- and c-quark decay in the Sun and Earth and include the effects of 
hadronization and stopping of heavy quarks and stopping and absorption of 
neutrinos in the Sun. They therefore can be used in place of Monte Carlo 
calculations. 



4. Neutrinos From Gauge Bosons 

A W boson decays directly to a muon and a muon neutrino 10.5% of 
the time, and a Z decays to a muon neutrino-antineutrino pair 6.7% of the 
time. Additional muon neutrinos are indirectly produced by W decay to h 
and c quarks and r leptons. The energies of these neutrinos are generally 
smaller, so the contribution provided by these indirect neutrinos to (^Nz'j is a 
small fraction of that from the direct neutrinos. Recall that light-fermion final 
states produce only very low energy neutrinos since they become thermalized 
before decaying [3], and their contributions to the neutrino energy moments are 
completely negligible. 

We will first calculate the i^N z'^j from decay of polarized and unpolarized 
gauge bosons in the Earth, and show that polarization is never more than a 
10% effect. Given these results, we will present expressions for dN/dE from 
decay of unpolarized gauge bosons only in the Sun and in the Earth, and then 
we will provide an analytic result for (^Nz'j for gauge-boson decay in the Sun. 

In the rest frame of the vector boson, the neutrino is emitted with an 
energy equal to half the vector-boson mass; in the laboratory frame, the energy 
of the emitted neutrino is Ei, = £"2(1 -|- /3cos6')/2 where 6 is the angle in the 
rest frame between the neutrino direction and the direction of motion of the 
center of mass, Ei is the energy at which the vector boson is injected, and /3 is 
its velocity. Therefore, for a vector boson injected into the Earth with energy 
Ei, the second moment of the resulting neutrino spectrum is 
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(4.1) 



where P(cos6) is the decay angular distribution in the rest frame of the gauge 
boson, and and F^ is the fraction of vector-boson decays which produce a neu- 
trino. 

For unpolarized gauge bosons, the decay is isotropic, and [5] 
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If the gauge boson has helicity /i = ±1 (transversely polarized), then 
P(cos9) (X (1 ± cos 6*)^, and 
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(4.3) 



and if the gauge boson is longitudinally polarized (/i = 0), then P(cos6) ex sin 6 
and 

{Elh = f (l + y) • (4-4) 

The velocities of WIMPs in the Sun or halo are t; ^ 1, so annihilation pro- 
ceeds only through the lowest angular-momentum state, the s wave. Although 
the s-wave cross sections for annihilation of neutralinos into ti [19][20][21], 
gauge-boson pairs [10], and Higgs-boson-gauge-boson pairs [5] were calculated 
previously, the helicity-amplitude formalism used by Drees and Nojiri [13] is 
needed to describe the polarization state of the annihilation products. Accord- 
ing to the results of [13], the gauge bosons produced in the s-wave annihilation 
processes XX ~^ W'^W~ and XX ~^ ZZ are produced with a tensor polariza- 
tion; that is, the gauge bosons are transversely polarized, half with spin aligned 
(/} = 1) with the direction of motion and half with spins antialigned (/i = — 1). 
Therefore, the scaled second moment of the neutrino spectrum from neutralino 
annihilation into W'^W~ pairs in the Earth is an average of that from the 
two transverse-polarization states. The second moment for neutrinos from the 
Earth is then given by 
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(4.5) 



where Vw -^ nv = 0.105. There will be additional muon neutrinos from r 
leptons from W decay, but their contribution to i^N z'^j is small. Only one 
in five r decays produces a muon neutrino, and the energy of this neutrino is 
typically only one third the energy of a neutrino from direct decay. Therefore, 
Eq. (4.5) properly accounts for polarization and should be accurate to better 
than 3%. 

The expression for ZZ pairs is obtained similarly to the above. We include 
an overall factor of two, which counts the Z bosons, and arrive at the result 



[Nz 



2\® 



1 



ixx^zz 



2Fz^.^, -(1 + -/32), 



(4.6) 



where Vz^v v = 0.067. Again, there will be additional neutrinos from Z 
decay to hh, and to a lesser extent, cc and fr, and the magnitude of their 
contribution to i^N z'^j is easily estimated. Although the branching ratio to hh 
is 2.2 times that into VfiVp,, a muon neutrino is produced in only one in ten 
&-quark decays. In addition, the &-quark energy is reduced to 0.7 times its 
original energy during hadronization, and the resulting muon neutrino carries 
only a third of that energy. By performing similar estimates for the r lepton 
and c quark, we find that Eq. (4.6) is accurate to better than 4%. 

Gauge bosons that are produced in the annihilation processes xX ~^ ZH , 
where H is one of the two neutral scalar Higgs bosons, or xX ~^ W'^ H~ (or 
its charge conjugate), where H is the charged Higgs boson, are longitudinally 
polarized [13], so Eq. (4.4) is the appropriate result to use for these final states. 
Thus we have 

(^•^')®x-Zff = Tz^.,.-, \{l + i/3|) + Higgs decay contribution. (4.7) 

The Higgs-decay contribution can be evaluated using the results of Section 2. 
Generally, the Higgs contribution is much smaller than that from the Z; there 
is an additional step in the decay chain, so the energies of the neutrinos are 
smaller. This is especially true if the Higgs is light enough that it decays only 
to light fermions. 

Previously, it was assumed that the gauge bosons were unpolarized, as 
in Eq. (4.2) [5]. Given our new results, we see that a proper treatment of 
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polarization changes the second moments by no more than 5% for the tensor 
polarization, and no more than 10% for the longitudinal polarization. In the 
rest of this Section, we will ignore the small effect of polarization and treat 
decays of unpolarized gauge bosons only. 

The neutrino energy distribution from W or Z bosons at rest is simply a 
delta function with support at half the gauge-boson mass. 
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where ms is the gauge-boson mass and F is the branching ratio for direct decay 
into muon neutrinos (this expression should be multiplied by 2 for Z bosons). 
Applying Eq. (2.1), the moving-frame distribution from unpolarized W or Z 
bosons moving with energy Ei and velocity /3 is 

(^^.^^^)^|r(mB7/3)-i forf (l-/3)<i?, <f (1 + /3), 
i otherwise. 

(4.9) 

Here we have neglected the contribution from other W decay modes. Although 

the contribution to dN/dE from these other modes may be significant at low 

energies, these low-energy neutrinos provide little contribution to the upward- 

muon fiux, as discussed above, so the final result obtained using Eq. (4.9) should 

be accurate to better than 4%. 

The energy distribution of neutrinos from gauge-boson decay in the Sun 

can be obtained by applying Eq. (2.6). The result for the scaled second moment 

of the neutrino distribution from W decay in the Sun is then [5] 

E=E,(l+l3)/2 
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(4.10) 
for VF's injected into the core of the Sun with energy Ei and velocity /3. 
Once again, the expression for ZZ pairs is obtained by replacing Tw^fiv with 
2Tz^i, V , where the factor of two counts the Z bosons. This does not include 
the contributions from the decays of the W and Z bosons to heavy fermions, 
but again, we have already seen that this effect is less than about 4%. 



Proper treatment of the effects of vector-boson polarization on the scaled 
second moment from the Sun leads to a very small correction to Eq. (4.10). In 
the low-energy limit, the effect of the polarization vanishes, as we have seen. 
In addition, the effect of the polarization must vanish in the high-energy limit 
as well. This is because the Lorentz boost factor suppresses the contribution 
of neutrinos with cos 61 ~ — 1 (antialigned with the gauge-boson direction of 
motion), but the interactions suppress those with cos 61 ~ 1 (aligned) because 
of their higher energy. The net result is that the integrated effect of the po- 
larization is never more than 2%, and is thus negligible for neutrinos from 
gauge-boson decay in the Sun. 

5. Neutrinos from Top-Quark Decays 

In many cases where the WIMP is heavier than the top quark, WIMPs 
annihilate predominantly to top quarks, especially if the top is as heavy as 
recent evidence suggests [22] [23]. Therefore, the neutrino spectrum from top- 
quark decay is needed for a proper calculation of the fiux of energetic neutrinos 
from WIMP annihilation in the Sun or Earth. 

In previous work [3] [5], the top-quark mass was assumed to be 60 GeV. In 
this case, calculation of the hadronization and decay channels of the top were 
quite difficult, and the resulting neutrino spectrum, calculated using Monte 
Carlo techniques, was accordingly uncertain. In addition, the results were 
highly dependent on the assumed top-quark mass. Given that the top quark 
is significantly heavier than the W [22] [23], the calculation must be repeated. 
The top quark is heavy, so it decays almost exclusively to a VF boson and a h 
quark. Hadronization effects are unimportant, and decay is essentially a free- 
particle weak decay. Therefore, our new results are less subject to theoretical 
uncertainty. 

The neutrino spectrum from top quarks follows from the neutrino spectra 
of its decay products, the W boson and the h quark, which we have already 
computed. We will first calculate the neutrino spectra from top-quark decay 
in the Earth and Sun and assume that the W bosons from top-quark decay 
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are unpolarized. We will then calculate (^Nz^j from top-quark decay in the 
Earth both with and without polarization and show that polarization of the 
gauge boson is only a 3% effect. Also, the results of Ref. [13] show that the 
positive and negative helicity amplitudes for top-quark production are equal. 
Therefore, the top quarks produced by neutralino annihilation are unpolarized, 
and we will consider unpolarized top quarks only. We will then present an 
analytic approximation to i^N z'^j for top-quark decay in the Sun. 

First consider the case where the top quark decays at rest. The top decays 
t -^ Wh with a branching ratio close to unity. The spectrum of neutrinos from 
top-quark decay at rest is given by a sum of the VF-decay contribution and the 
&-decay contribution: 
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In the rest frame of the top quark, the W has energy 
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and its velocity is 
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The VF-boson contribution is given by Eq. (4.9) with injection energy Ew ■ The 
h quark is first slowed from its injection energy, Ei, = (m^ — m'^)/(2mt), by 
hadronization to an energy Ed = ZfEi,, and has a velocity /3 close to unity. The 
&-quark contribution is then given by the (3^1 limit of Eq. (3.8). 

The boost formula, Eq. (2.1), yields the neutrino spectrum from unpolar- 
ized top quarks moving with velocity j3t and 7^ = (1 — j3'^)~^l'^ , 
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where e± = Ew{^ ± /3vi/)/2 with Ew and (3w equal to their values in the top- 
quark rest frame [Eqs. (5.2) and (5.3)], E± = E„j^ (1^ /3t)~^ , and Q(x) = 1 if 
X is true and Q(x) = otherwise. Also, the &-quark decay energy is Ed = ZfEi, 
where Ei, is the energy of the h quark in the top-quark rest frame. The function 
Di, is given by 



Di[x, y] = ]- \9(x^ - j/2) + 4(x^ - j/^) + 6 In ^ 
3 L X 



(5.5) 



The kinematics of the decay dictate that the h quarks carry a smaller fraction of 
the energy than the W bosons, their energy is further degraded by hadroniza- 
tion, and they undergo three-body — rather than two-body — decay. Therefore, 
even though the &-quark contribution to (dN/dEy^f*' may be significant at low 
energies, it will contribute only about 9% to the final result for the upward- 
muon fiux. 

The distribution from the Sun is obtained by applying Eq. (2.6) to Eq. 
(5.4). It should be noted that application of Eq. (2.6) to Eq. (5.4) does not 
take into account the stopping of h quarks in the Sun; however, neutrinos from 
h quarks contribute only 9% of the upward-muon signal, so little accuracy will 
be lost by neglecting this effect. 

We now calculate i^N z'^j for top-quark decay in the Earth, and quantify the 
effect of gauge-boson polarization on the results for the upward-muon signal. 
Consider the case where the top quark decays at rest. A fraction [24] 
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(5.6) 



of the W bosons are produced in the the longitudinal-helicity state, and the 
rest (a fraction 1 — f^) are produced in the transverse-helicity states, with 
equal probability for positive and negative helicity states. The result for (^Nz^j 
for top quarks at rest (Ei = nii where Ei is the top-quark energy) then follows 
from the results of the previous Section. The only subtlety is that the second 
moment must be scaled by the square of the top-quark mass rather than the 



21 



22 



square of the VF-boson energy. The contribution from the semileptonic h decays 
can also be included. Thus, we have 
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(5.7) 
For mt = m GeV [22], /i ^ 0.7 and /3^ = 0.4. If the W bosons were 
produced unpolarized, then we would have f^ = 0.33, so that consideration of 
the polarization of the W bosons in the decay W ^ jiv results in a small change 
from the unpolarized case, less than about 3%. And as discussed above, the 
contribution to (Nz'^\ from h quarks is subdominant (about 9% of the total). 

The top quark decays immediately, so the scaled second moment from 
decay of a top quark in the Earth moving with velocity (3t and energy Et is 
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Therefore, the i^N z'^j for tops moving with velocity j3t is just (1 + 13^/3) times 
that for top quarks at rest, as given in Eq. (5.7). 

Eqs. (5.7) and (5.8) provide expressions for the scaled second moment 
of the neutrino distribution from top quarks from WIMP annihilation in the 
Earth, accurate to at least 3%. 

In order to obtain (Nz'^\ from top-quark decay in the Sun, the stopping 
and absorption of neutrinos, as given in Eqs. (2.4) and (2.5), and stopping of 
the &-quark hadrons must be included. 

Neutrinos come predominantly from the W boson produced by the top- 
quark decay, and to a lesser extent from the h quarks. We already know the 
i^N z'^j from injection of these particles in the center of the Sun. The only catch 
is that, to use Eq. (4.10), we must ignore polarization of the W; the results of 
the previous subsection show that W polarization in top decay is less than a 
3% effect. Therefore, the easiest way to obtain the result for (Nz"^) from top 



quarks from the Sun is to integrate the (^Nz^j for W bosons and b quarks over 
the appropriate injection energies. 

The expression for (Nz'^\ as a function of the top-quark energy, Et, for 
top quarks injected into the Sun follows from Eq. (2.12) and is 
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where we have included the sum over both the W bosons and h quarks from top- 
quark decay. The integral is over the injection energy of the decay particles, and 
the moments {Nz^)f include the effects of interactions with the solar medium. 
For W bosons, i^N z'^j ' is given by Eq. (4.10), and for h quarks by Eq. (3.15). 
The quantities Et, (3t, and ■jt are the energy, velocity, and 7 factors of the 
top quark, and Ew and (3w are the energy and velocity of the W boson in 
the rest frame of the top, as given in Eqs. (5.2) and (5.3). The quantities 
Ell = (toj — m'^)/(2mt) and /3j ~ 1 are the energy and velocity of the h quark 
in the rest frame of the top. Also, recall that i^N z'^j for neutrinos from the Sun 
are different than those for antineutrinos. 

The results for (Nz'^\ for neutrinos and antineutrinos from injection of 
top quarks into the core of the Sun and Earth are shown in Fig. 1. There are 
three bands of curves: the top is the result for (^Nz^j from the Earth (and 
is the same for both neutrinos and antineutrinos). The center band is the 
result for antineutrinos from top quarks injected into the core of the Sun, and 
the lower band is that for neutrinos from top quarks injected into the core of 
the Sun. The curves for antineutrinos from the Sun are higher than those for 
neutrinos since neutrinos interact more strongly with the solar medium than 
antineutrinos and are therefore more strongly attenuated as they pass through 
the Sun. Comparison of the results for (Nz'^\ from the Sun and Earth illustrates 
that absorption and stopping of neutrinos in the Sun are very important effects. 

In each band, the solid curve is the total result for i^N z'^j which accounts for 
neutrinos from both the W boson and h quark. In each case we have taken nit = 
174 GeV [22]. The dashed curve is the result obtained ignoring the contribution 
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Fig. 1. Plot of (Nz'^\ for top quarks injected into tlie core of tlie Sun and 
Eartli as a function of tlie injection energy Et- See text for a description of the 
curves. 

from h quarks. For neutrinos from the Sun, the &-quark contribution is about 
10% at low energies, and increases to roughly 25% at energies around 1000 
GeV. The increased importance of the &-quark neutrinos from the Sun at higher 
energies is due to the fact that the neutrinos from W bosons are more energetic 
and therefore more strongly absorbed as they travel through the Sun. 

In Fig. 1, we also illustrate the effect of varying the top-quark mass. The 
dotted curves above and below the solid curves in each band are the total (N z"^) 



obtained using nit = 131 GeV (a conservative lower limit [23]) and nit = 225 
GeV (the 3(T upper limit suggested by the CDF data [22]), respectively. 

We have crafted the following analytic fit which provides a good approxi- 
mation (better than 10%) to the (Nz'^\ over the range nit < Et < 3000 GeV 
for neutrinos from the Sun (the two lower solid curves curves in Fig. 1): 
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For neutrinos, Ai, = —0.825, 5^ = —3.31, and C'l, = —5.39. For antineutrinos, 
Ap = —0.889, Bp = —3.94, and Cp = —6.40. Of course, if the top-quark mass 
differs from 174 GeV, these coefficients will change slightly. 

6. Summary 

In this paper we have calculated the differential energy spectra of muon 
neutrinos and antineutrinos from particle decay in the Sun and Earth. These 
results will be useful for computing the rates for upward-muon events in neu- 
trino telescopes from annihilation of WIMPs that have accumulated in the Sun 
and/or Earth. The analytic expressions we have derived for neutrino spectra 
from decay of particles injected in the Sun result from the two- or three-body 
decay kinematics and include the effects of hadronization. The results for par- 
ticle decay in the Sun take into account the additional effects of stopping of 
heavy hadrons in the Sun and stopping and absorption of neutrinos as they 
pass through the Sun. Neutrino spectra are given for all the WIMP-annihilation 
channels that will produce energetic neutrinos. These expressions will give re- 
sults for the fiux of neutrino-induced upward muons accurate to better than 
0(10%) and can be used in place of Monte Carlo calculations. 

In particular, the energy spectra from injection of tt, cc, and hh, pairs in 
the Earth are given by Eqs. (3.2), (3.8), and (3.11), respectively. The neutrino 
spectra from injection of W'^W~ and ZZ pairs in the Earth are given by Eq. 
(4.9), and the spectrum from injection of top quarks in the Earth is given by 
Eq. (5.4). The antineutrino spectra from particle decays in the Earth are 
the same as the neutrino spectra. The neutrino spectra from injection of rf, 
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W'^W~ , and ti pairs in the Sun are obtained by applying Eq. (2.6) to Eqs. 
(3.2), (4.9), and (5.4), respectively, and the neutrino spectra from injection 
of hh and cc pairs in the Sun are given by Eqs. (3.9) and (3.12), respectively. 
Recall that interactions of antineutrinos with the solar medium differ than those 
of neutrinos, so the antineutrino and neutrino spectra from particle decay in 
the Sun differ. The neutrino spectra from decays of Higgs bosons in the Sun 
and Earth can be obtained from the spectra of the Higgs decays products by 
applying Eq. (2.11). 

The probability of detecting a neutrino-induced upward muon is propor- 
tional to the square of the neutrino energy; one power comes from the charged- 
current cross section for muon production, and the other comes from the muon 
range. We argued that in most cases, the detector thresholds are small enough 
compared to the WIMP mass that the thresholds can be ignored. If so, then 
the upward-muon flux from WIMP annihilation is proportional to the second 
moment of the neutrino energy distribution, Eq. (1.2) will provide an accurate 
estimate of the muon flux from a given WIMP candidate, and detailed informa- 
tion about the shape of the neutrino energy distribution is not required. For this 
reason, we have provided analytic expressions for the scaled second moments 
of the neutrino distribution, (#z^), for each WIMP-annihilation channel. The 
scaled second moments from injection of rf, hh, cc, W'^W~ , ZZ, and ti pairs 
in the Earth are given by Eqs. (3.3), (3.13), (3.14), (4.5), (4.6), and (5.8), 
respectively. The scaled second moments from injection of rf and ti pairs in 
the Sun are given by Eqs. (3.4) and (5.10), respectively. Those from injection 
of hh and cc pairs in the Sun are given by Eq. (3.15), and those from injection 
of W'^W~ and ZZ pairs are given by (4.10). The scaled second moments from 
Higgs-boson decay in the Sun and Earth can be obtained from those of the 
Higgs decay products using Eq. (2.12). 



The results for the 



id moments 



, ,,6^ 



(Nz^), for neutrinos and an- 



tineutrinos from injection of particle-antiparticle pairs in the Sun with energy 
equal to the WIMP mass, m^, are plotted in Fig. 2 and Fig. 3, respectively. 
The strongest signals are from those particles that decay immediately and di- 
rectly to neutrinos, the tt and gauge-boson final states, and the signal from top 
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Fig. 2. The second moment, m^^ (Nz'^\, of the neutrino energy distribution 
from injection of particles with energy equal to the WIMP mass m^ in the Sun. 
The solid curve is for tt pairs, the upper (lower) dashed curve is for W'^W~ 
{ZZ) pairs, and the upper dot-dash curve is for rf pairs. At the bottom are 
the hh and cc curves, and the hh curve is slightly higher than the cc curve. 

quarks is only slightly weaker. The turnover in the gauge-boson signal strength 
is due to absorption of high-energy neutrinos in the Sun; the signals from other 
particles will turn over similarly at larger energies. The Figures illustrates the 
importance of hadronization and stopping of heavy hadrons: the signals from 
h and c quarks are significantly smaller than those from particles that do not 
from hadrons. 

The results listed here contain the effects of hadronization of heavy quarks 
and interactions of heavy hadrons and neutrinos with the solar medium. We 
have included only the neutrinos from prompt particle decay. In some cases, 
there will be additional neutrinos from secondary decays, but their contribution 
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Fig. 3. Same as Fig. 2, but for antineutrinos. 

to the upward-muon flux will be negligible. For example, there will be neutrinos 
produced by gauge-boson decay to r leptons and h and c quarks in addition 
to the prompt neutrinos from gauge-boson decay, but we estimate that their 
contribution to the upward-muon flux will be less than 4% of the total. Also, 
there will be additional neutrinos from the c quarks produced by h decay, but 
we estimate that these will contribute roughly 10% of the total upward-muon 
flux. The additional contribution to the upward-muon flux from decays of r 
leptons produced by charged-current interaction of r neutrinos as they pass 
through the Sun are negligible. If our neutrino spectra are used to compute 
event rates for detection techniques which are proportional to the first moment 
of the neutrino energy distribution (such as searches for contained events), then 
the contribution from the neutrinos from secondary decays which we neglect 
may be significant and the event rates will be underestimated. 



The neutrino energy spectrum from decay of a given particle depends on 
the energy of the decaying particle, and it may depend on the polarization state 
as well. We have evaluated the effect of polarization on the upward-muon fiux; 
we found that it is never more than a 10% effect and can typically be accounted 
for easily 

The largest theoretical uncertainty involves modeling the effect of stopping 
of heavy hadrons in the Sun. This can be traced to the uncertainty in the values 
adopted for the stopping coefficient, E^ [c.f. Eq. (2.7)]. However, neutralino 
annihilation into h and c quarks is always accompanied by annihilation into 
r leptons (with a branching ratio enhanced by radiative corrections [25]). If 
heavy enough, WIMP annihilation is predominantly to ti pairs and possibly 
gauge-boson final states, so the uncertainty in the total upward-muon fiux due 
to uncertainty in the stopping effect is diluted. For all annihilation channels, 
there is also the uncertainty in the absorption and stopping length for neutrinos 
and antineutrinos in the Sun, but this is small at energies near a TeV, and 
is virtually negligible at lower energies. These uncertainties do not effect the 
neutrino spectra from the Earth. Note that these theoretical uncertainties must 
also enter Monte Carlo calculations and are not introduced by our analytic 
treatment. 
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